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Abstract 

We study the influence of unitarity corrections on the Drell-Yan transverse momen- 
tum distribution within the color dipole approach. These unitarity corrections are 
implemented through the multiple scattering Glauber-Mueller approach, which is 
contrasted with a phenomenological saturation model. The process is analyzed for 
the center of mass energies of the Relativistic Heavy Ion Collider (RHIC, y/s = 500 
GeV) and of the Large Hadron Collider (LHC, yfs = 14 TeV). In addition, the re- 
sults are extrapolated down to current energies of proton-proton collisions, where 
non-asymptotic corrections to the dipole approach are needed. It is also shown that 
in the absence of saturation, the dipole approach can be related to the QCD Compton 
process. 

PACS numbers: 13.85. Qk; 12.38.Bx; 12.38.Aw. 



1 E-mail : mandr ebc@if . ufrgs . br 
2 E-mail:gay@if. ufrgs.br 
3 E-mail:magnus@if. ufrgs.br 
4 E-mail:jorgr@lanl.gov 



2 



I. INTRODUCTION 

The high energies available in the hadronic reactions at RHIC (BNL Relativistic Heavy 
Ion Collider) and to be reached at LHC (CERN Large Hadron Collider) will provide a 
better knowledge concerning parton saturation. In such a kinematical region the production 
of massive lepton pairs in hadronic collisions (Drell-Yan (DY) process |l|) can be used to 
investigate the high parton density limit, since it is a clean reaction probing the gluon 
distribution through the QCD Compton process. In particular, the Drell-Yan transverse 
momentum (j>t) distribution can be expected to be sensitive to saturation effects. 

Saturation and nuclear effects are most conveniently described within the color dipole 
approach which is, in fact, especially suitable for this purpose (see Refs. 0, U, O, E| 
for some applications). The dipole approach is applicable only at high energies, and it is 
formulated in the target rest frame, where the DY process looks like a bremsstrahlung of 
a virtual photon decaying into a lepton pair (see Fig. [TJ. The advantage of this formal- 
ism is that the DY cross section can be written in terms of the same color dipole cross 
section as small-x Deep Inelastic Scattering (DIS). Although diagrammatically no dipole is 
present in bremsstrahlung, the dipole cross section arises from the interference of the two 
bremsstrahlung diagrams, see Ref. |7| for a detailed derivation. The cross section for a ra- 
diation of a virtual photon from a quark scattering on a nucleon (N) can be written in a 
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FIG. 1: In the target rest frame, DY dilepton production looks like a bremsstrahlung. A quark or 
an anti-quark from the projectile hadron scatters off the target color field (denoted by a curly line) 
and radiates a photon (7*) with mass M (before or after the quark scatters), which subsequently 
decays into the lepton pair 



factorized form as 

da T)L (qN -> 7*X) 
din a 



J d 2 r ± \^ 1 y^(a,r ± )\ 2 a d ip(ar ± ), (1) 



where <idi V is the same dipole cross section as in DIS, which should take into account non- 
perturbative and saturation effects at high energy |3j. The energy dependence of a^ip is 
not explicitly written out. Here, r± is the photon-quark transverse separation, and the 
argument of the dipole cross section, ar±, is the displacement of the projectile quark in 
impact parameter space due to the radiation of the virtual photon, different to the DIS 
case, where the dipole separation is just r±. The ^j* q are the light-cone wave functions 
for radiation of a transversely (T) or longitudinally (L) polarized photon (see e.g. Ref. []| 
for explicit expressions). While the light-cone wave functions are calculable in perturbation 
theory, the dipole cross section can be determined only with input from experimental data. 

The goal of this work is to investigate the influence of unitarity corrections on the DY 
dilepton px distribution, describing these unitarity corrections by the multiple scattering 
Glauber-Mueller approach [8| and including them into the dipole cross section. The results 
are contrasted with the QCD improved phenomenological saturation model of (Tdipi Ref. El> 
which quite successfully describes DIS and diffractive DIS data. In lines of a previous work 
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here it is investigated the role of the 7*g wave functions in the px distribution, charac- 
terizing the relation between dipole sizes and transverse momentum. A striking advantage 
of the color dipole picture is a finite cross section for the lepton pair distribution at small 
Pt — > 0, even in the leading order calculation, feature associated with the saturation encoded 
in the dipole cross section. In the conventional parton model, the perturbative calculation 
of 0(a s ) yields a divergence at pt = 0, and one has to resume large logarithms, In (p T /M 2 ), 
in an appropriated scheme fiol ]. in order to obtain a physically sensible result. 

The large amount of low Bjorken-a; DIS available data allows to constrain the dipole 
cross section at very high energies and to calculate the DY cross section without additional 
free parameters. However, in the current energies of the hadronic colliders there are non- 
asymptotic corrections to the dipole cross section, which have to be taken into account, in 
order to describe experimental data. Therefore, one also introduces a parameterization for 
that contribution, which is negligible already at RHIC energies. In addition, we show how 
the dipole approach for the DY px distribution is related to the QCD Compton process, 
which contributes at order a s to the conventional parton model of DY dilepton production. 



FIG. 2: Production of massive photons through the QCD Compton process. The subsequent 
decay of the 7* (wavy line) into the dilepton pair is not shown here. Curly lines denote gluons, 
quarks are represented by lines with arrows. 

The two approaches are equivalent in a certain limit. 



II. RELATING DIPOLE APPROACH AND PARTON MODEL OF HIGH p T 

DILEPTON PRODUCTION 



Although the dipole approach and the next-to-leading order (NLO) parton model have 
been compared numerically in Ref. JjjJ, one may still wonder, how these two approaches can 
be related to one another analytically. This will be the topic of the present section. Using 



the leading order expression 



III, 



7r 2 a. 



r\ x G(x) 



(2) 



for the dipole cross section, it can be demonstrated how the dipole approach for high p? 
dilepton production is related to the QCD Compton process, see Fig. El In Eq. (0), xG{x) is 
the density of gluons with momentum fraction x in a nucleon, and a s is the strong coupling 
constant. First, we shortly review the formulas for high dilepton production in the dipole 
approach and in the parton model, before we show how they can be translated into each 
other. 

In the dipole approach, the DY transverse momentum distribution is given by [12], 
d 3 a(pp -> 1+rX) 



dydM 2 dp 2 T 



3M 2 



X\ 



da 



Oi 



Q 



x 



q 



x x 



a J \ a 

' 'a, r ± )C* (a, r' ± ) + r ± )*%Ja, r' ± ) 



- 7«g 



x - [<J dip (x, ar±) + a dip (x, ar' ± ) - a dip (x, a \r± - 



(3) 
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where the quark (antiquark) distributions in the projectile are denoted by q (q respectively). 

The usual definitions of the kinematic variables are employed, i.e., 

2iVg 2P 1 -q 

xi = , x 2 = , (4) 

s s 

where q is the four momentum of the virtual photon (M 2 = g 2 ), and P\ 2 are the four 
momenta of the projectile (1) and target (2) hadron. By evaluating the scalar product for 
X\ in the target rest frame, it is easy to show, that the projectile parton distributions in 
Eq. (J2J) are probed at momentum fraction Xi/a, where a is the momentum fraction taken 
by the photon from the projectile quark. Furthermore, Pt is the transverse momentum of 
the 7* in a frame with the z-axis parallel to the projectile quark, and 

-Ml) < 5 > 

is the rapidity of the photon. In addition, 

rj 2 = (1 - a) M 2 + a 2 m 2 g . (6) 

The quark mass m q is set to in this section. The upper limit of the a-integration in Eq. (jSJ) 
is determined from the condition that the invariant mass of the final state cannot exceed the 
total available center of mass (cm.) energy of the projectile quark-target nucleon system, 
i.e., 

xis p 2 T + rf _ _ p 2 T 

a a (1 — a) x\s — M l 

where \fs is the hadronic cm. energy. In the high energy approximation, a max = 1 for 

s — > oo. In this section, however, we work with the exact value of a max . 

With <j dip given by Eq. (0), the integrals over r± and r' ± in Eq. (jSJ) can be performed 

analytically with the result |l2j, 

fd 3 a(pp -> l+l~X)\ 



dydM 2 dp\ 



-approx. 



2 ^max N f 



< ^^x l I daY, e 2 a 
9M 2 J ^ q 



xG(x) 



1 + (1 _ a) 2 } A±±_ + 4M 2 (1 _ a f PL_ I . (8) 

v J [p 2 T + r] 2 ) (Pt + V 2 ) J 

In order to obtain Eq. (jSJ), one has to assume that xG(x) does not depend on r±_ through 
scaling violations. Note also, that the r^-approximation, Eq. (J2J), is applicable only at large 
Pt- 
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In the parton model, on the other hand, the high distribution of DY dileptons produced 
via the QCD Compton process, see Fig. El is given by 



f d 3 a(pp^l + l-X) \ 
\ dydM 2 dp\ J 



Compton 



2 

9M 2 



x 



d,Xn 



t 



X a Xi) 



N f 



e l {[<? ( x a) + 1 ( x a)] G(x b ) + G(x a ) [q (x b ) + q (x b )}} 



q=l 



u 



-2M Z — - - - - 

su u 



(9) 



(see e.g. Ref. [13j for details). In Eq. @, x a and x b are the momentum fractions of the 
colliding partons, and 



x„ 



x 1 - M 2 /s 



x 2 



(10) 



Note that at finite p?, x ajb ^ Xi j2 , where Xi >2 = y T+ s M e^ +, ^ y . The partonic Mandelstam 
variables, s, i, u, are defined in terms of x a , x b and the four-momenta of the colliding 
hadrons, see Fig. |2j In order to compare Eqs. (jHJ) and ®, one has to express the partonic 
Mandelstam variables in terms of a and p 2 ^, 



Furthermore, 



s = {x a Pi + x b P 2 ) 
i=(q- x b P 2 f 
u = {q- X a P{f 



- £l 

X n , Xjy 



Pt + V 2 
a (1 — a) 

Pt 



1 — a 

p T + r] 2 



a 



p T + r] 2 



x 2 . 



(11) 
(12) 
(13) 

(14) 



a (1 — at) p\ + rf 

Inserting the expressions for the partonic Mandelstam variables, Eqs. (JHJ), (fl2|) and (jH 
into Eq. (jHJ), one obtains a result very similar to Eq. (JHJ), except for the combinations of 
parton distributions, 



f d 3 a(pp^l + l-X) \ 
{ dydM 2 dp\ J 



9M 2 



-Xi 



■t-i 



Compton 

da^2e 2 q {[q (x a ) + q (x a )] x b G(x b ) + G(x a )x b [q{x b ) + q (x b )}} 

9=1 



X 



Pt + V +4M 2 (l-«) 2 



Pt 



(Pt + V 2 ) 



(Pt + V 2 ) 



(15) 
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When saturation effects are neglected, the dipole approach reproduces that part of the 
QCD Compton contribution to DY, in which the quark comes from the projectile and the 
gluon from the target. Thus, the dipole approach is valid, when the first term in the convo- 
lution of parton distributions in Eq. ()15j) dominates. This is the case at large rapidity and at 
small Xb, both conditions fulfilled. The range of validity of the dipole approach can of course 
only be established a posteriori. This is similar to the problem of determining the lowest 
scale at which perturbative QCD still works. The dipole approach is phenomenologically 
successful for values of x 2 < 0.1, though most parameterizations of the dipole cross section 
are fitted only to DIS data with Bjorken-x < 0.01. 

Regarding the rapidity (y) range, in which the dipole formulation can be applied, some 
guidance on the minimal value of y can be obtained from the numerical comparison of the 
dipole approach and the next-to-leading order (NLO) parton model in Ref. |7J. At RHIC 
energy ^/s = 500 GeV, virtually no deviations between the dipole approach and the NLO 
parton model have been found for y > 0.5 [7]. This means that one can safely compare the 
dipole approach to future DY measurements from the two PHENIX muon arms . 

On the other hand, the dipole approach takes into account several effects that will be 
important at high energies. A realistic parameterization of the dipole cross section includes 
gluon saturation, which is not contained in the standard parton model. Moreover, a dip 
contains information about the transverse momentum distribution of the target gluons, 
thereby is more complete than the gluon distribution in the collinear factorization approach. 
Finally, with a realistic parameterization of the dipole cross section at large separation r±, 
one can apply Eq. (JHJ) also at low pr, while the conventional parton model pr distribution, 
Eq. (JHJ), is applicable only at very high p T > M. 

III. FEATURES OF THE DRELL-YAN CROSS SECTION IN THE COLOR 

DIPOLE APPROACH 

In this section, it is investigated which distances r_|_ in impact parameter space are impor- 
tant for the Drell-Yan pt distribution. For this purpose, the behavior of the weight function 
for a di V as function of p = ar± for different values of pt is studied. 

Three of the four Fourier integrals in Eq. (JHJ) can be performed analytically with the 
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result 



a, 



da 



DY 



1 



dM 2 dx F d 2 p T 6 vr 3 M 2 (x x + x 2 ) -A) 
where the weight function W(p, pr) is given by 

" l da x\ A 



dpW(p,p T )a dip (p), 



(16) 



x <j [m^a 4 + 2M 2 (l -a) 2 ] 

VPt 



+ [1 + (1 - a) 



— ^ t i(p) - -r T 2{p) 
Pt + 7 4ry 

rp ( \ T l(P) ,Vrrf\ 

2 | 2 T 3 P ^ + T T 2 P 



(17) 



and the functions T,- read, 



Ti(p) = pJ (p T p/a)K (r]p/a)/a, 
T 2 (p) = p 2 J (pTp/a)K 1 (r]p/a)/a 2 
T 3(p) = pJi(pTP/a)Ki(rjp/a)/a. 



(18) 
(19) 
(20) 



The functions Jo and J\ are the first class Bessel functions of order and 1, whereas K Q and 
K\ are the second class modified Bessel functions of order and 1 (MacDonald functions). 

It was shown in Ref. js] that for the (p? integrated) mass distribution, the wave functions 
select the small p region. Large values of p > 2/M are exponentiated by the functions K 01 . 
It should be stressed here that large dipole sizes correspond to the non-perturbative sector 
of the reactions, whereas small size configurations give the perturbative piece. 

In the particular case of the dilepton p? distribution, a different picture is designed. In 
Fig. 01 we show W(p,pt) as a function of the photon-quark transverse separation p for 
typical fixed lepton pair mass M = 6.5 GeV and Feynman-x {xp = 0.625). The results are 
presented for two center of mass energies: the plot on the left corresponds to \fs = 38.8 
GeV (available at the E772), whereas in the plot on the right, y/s = 500 GeV (RHIC). For 
the effective light quark masses, the value m q = 0.2 GeV was used. Three different values 
for the dilepton transverse momentum were selected, px = 0, 1 and 4 GeV. 

As can be seen from Eq. f!17j1 . the oscillating Bessel functions J, drive the behavior of 
the W(p,pt) as a function of p. The following general picture can be drawn from the plots: 
for large px the large dipole size configurations get suppressed, because W(p,px) is rapidly 
oscillating. This suppression mechanism is different from the exponential suppression of large 
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FIG. 3: The weight function W(p,px) as a function of p for different pj> at fixed xp = 0.625 and 
M = 6.5 GeV. 

dipole sizes in the case of the px integrated cross section, and complicates the numerical 
calculation of the pp distribution. On the other hand, as pp decreases, large p configurations 
become more important. The case pp — is of particular interest, since the weight function 
W(p,pt) selects very large dipole configurations and such a region is enhanced by increasing 
the energy. Therefore, the non-perturbative sector of the process should drive the small pp 



regime. On the other 
dipole configurations 



land, the large pt behavior is almost completely dominated by small 
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. These features are exploited in the next section, where are also 



discussed the different models that were employed for the dipole cross section. 



IV. THE DIPOLE CROSS SECTION 



The cross section for a small color dipole scattering on a nucleon can be obtained from per- 



turbative QCD . However, there are large uncertainties stemming from non-perturbative 
effects (infrared region) as well as from higher order and higher twist corrections. In the 
leading ln(l/x) approximation, the dipole interacts with the target through the exchange 
of a perturbative Balitsky-Fadin-Kuraev-Lipatov (BFKL) Pomeron, described in terms of 

n n 

the ladder diagrams [Iq . In the double logarithmic approximation, the BFKL equation |lfi| 
agrees with the evolution equation of Dokshitzer et al. (hereafter DGLAP equation). In 



10 



this limit, the dipole cross section reads, 

<r*p(x, r ± ) = ^rlx G DGLAP (x, Q 2 ) , (21) 

where xG DGLAP (x, Q 2 ) is the usual DGLAP gluon distribution at momentum fraction x and 
virtuality scale Q 2 = \/r\. The factor A appearing in the virtuality scale Q 2 = X/r\, has 
been taken a, A = 4 g althongh sam e .nagnitude va!ue S a re e q uivnlent at leading logarith- 
mic level 18 . The main feature of the dipole cross section above is the color transparency 
property, i.e., a dip ~ r\ as r±_ — ► 0. At large dipole size, the dipole cross section should 
match the confinement property a dip ~ &o- Concerning the large transverse separation (non- 
perturbative sector), our procedure is to freeze the r\ in Eq. (|2*Tj) at a suitable scale larger 
than r 2 nt , which corresponds to the initial scale on the gluon density perturbative evolution, 

Ql = 4/r c 2 ut - 

At high energies, an additional requirement should be met: the growth of the parton 
density (mostly gluons) has to be tamed, since an uncontrolled increasing would violate the 
Froissart-Martin bound, requiring the black disc limit of the target has to be reached at 
quite small Bjorken x. This feature can be implemented by using the multiple scattering 
Glauber-Mueller approach (GM), which reduces the growth of the gluon distribution by 
eikonalization in impact parameter space |8J. Therefore, one substitutes xG DGLAP in Eq. (JHJ) 
by the corrected distribution including unitarity effects, x G . A more extensive derivation 
of the GM dipole cross section and the expression of x G GRI can be found in the Sec. Ill 
of the Ref. [^. Following previous work Q], one shall use x = Xi as the energy scale in the 
dipole cross section, since x<i in DY is the analog of Bjorken-x in DIS. Note that x = axi 
was used in [7], however, the factor a has only a small numerical influence. 

Once the dipole cross section is known, one can also calculate the DY differential cross 
section, Eq. (fT7)|) integrated over p?, and compare it with the available data at small x-i- 
However, the current data on DY reactions are measured in a kinematical region where X2 
still takes rather large values, that is, x 2 — 0.1 at ^Js = 38.8 GeV, where the color dipole 
picture reaches the limit of its validity. Therefore, in order to compare the theory with these 
data, some procedure should be taken to extend the applicability of the dipole cross section 
at large x%. 

Note, that the dipole cross section Eq. (|2Tj) represents the asymptotic gluonic (Pomeron) 
contribution to the process, and at large x (low energy) a non-asymptotic quark-like content 
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should be included. In the Regge theory language, this means a Reggeon contribution, and 
therefore, we added the term 3J], 

a%=a rl X ^{l-xf. (22) 

to the dipole cross section, Eq. (J2H). 

Using the expre SS1 on above, good reSu .t S are obtained in describing the E772 data Q 
on mass distribution with a Reggeon overall normalization cq = 8 [3], reproducing similar 
results considering the saturation model Q|. Nevertheless, Eq. ()22j) has a shortcoming when 
one calculates the dilepton p T distribution: due to the fact that the weight function, Eq. (fTTj) . 
selects large dipole configurations at small pt (see discussion in the previous section) the 
~ r\ behavior in the Reggeon dipole cross section produces a non-negligible contribution 
at small pt even at RHIC energies. Therefore, Eq. ({22} was modified in order to cure this 
shortcoming and preserve our previous results. The Reggeon contribution now reads, 



fi v = <tqt]_x g va i (x, Q 2 ) , (23) 



where the quantity g va i is the valence quark distribution from the target and a reasonable 
description of the same E772 data is obtained with a value er = 7. The scaling violation 
from the valence parton distribution takes care of the steep growing on r±, which is present 
in the simple parameterization of Eq. (|22j) . removing the already mentioned shortcoming in 
the pt distribution at high energies. . 

Our main goal here is to investigate the DY p? distribution, using the GM dipole cross 
section. However, for sake of comparison, this analysis is contrasted with the phenomeno- 
logical saturation model of Bartels et al. (BGBK dipole cross section hereafter), Ref. jjj, 
which also includes the features of the dipole cross section discussed above. The model of 
Ref. is a QCD improved version of the saturation model of Ref. J^. The new model 
explicitly includes QCD evolution, and the dipole cross section is given by, 

Odi P {x, r±) = (7 1 1 - exp I — I > , (24) 

where the scale \j? is assumed to have the form 

V 2 = -T + f4- ( 25 ) 

n r± 

The authors of Ref. [£| propose the following gluon distribution at initial scale Ql = 1 GeV 2 , 

xg(x,Q 2 ) = A g x- x *(l-x) 5 - 6 . (26) 
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Altogether, there are five free parameters (a , C, /ig, A g and \ g ), which have been determined 
in Ref. [I| by fitting ZEUS, HI and E665 data with x < 0.01. In this fit the parameter o"o 
is fixed at 23 mb during the fits as in the original model, Ref. j^j. Here, we employ fit 1 of 
Ref. Q. 

In Ref. 0, where the old saturation model of Ref. 20] was used, the dipole approach 
was extrapolated to larger x<i by introducing a threshold factor into the saturation scale, 
i.e. Q 2 — » Q 2 (1 — X2) 5 . The factor (1 — x 2 ) 5 is motivated from QCD counting rules and 
suppresses the large x 2 contribution in the DY cross section. In our case, employing the GM 
or the BGBK dipole cross section, the large X2 threshold factor is already included in the 
collinear gluon distribution function. 

In addition, in Ref. Q|, (1 — x\)M 2 was used as the virtuality scale, at which the projectile 
parton distribution is probed, see Eq. (fTBj) . In this work, we shall use M 2 instead. The factor 
(1 — x\) is important only at large Xi, but has no effect at midrapidity. In the next section 
we study the dilepton transverse momentum distribution, making use of the results obtained 
above for low energies 



V. THE DILEPTON TRANSVERSE MOMENTUM DISTRIBUTION 

In this section, the DY dilepton transverse momentum distribution is calculated, using 
the Glauber- Mueller dipole cross section, Eq. (|2~T|) . and compared with the results obtained 
with the improved saturation model, Eq. (J2U). We will consider typical values for mass and 
xp. The projectile structure function employed was the LO GRV98 parametrization j^lj to 
the GM predictions and CTEQ5L (2^] for the saturation model ones. 

Before doing that, some comments are in order. The unitarity effects in the target will be 
significative at large rapidity y = 1/2 Ib.(xf/x2 + 1). In the central rapidity region (y ~ 0) 
the effects in the projectile could be also sizeable. In the last case, those effects in the quark 
distribution are smaller than in the gluon content. Therefore they will be disregarded in 
what follows. 

In Fig. |U results for the energies from RHIC (y/s = 500 GeV) and LHC (y/s = 14 
TeV) are shown with M = 6.5 GeV and xp = 0.625. At these energies and kinematics 
variables, the valence content is completely negligible. We emphasize that the xp value 
considered above, is an extreme case, where the rapidity variable acquires large values for 
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RHIC (y ~ 3) and LHC (y ~ 7) energies. In order to investigate the unitarity effects for this 
observable, the following comparisons are performed: The long-dashed curves are calculated 
with the dipole cross section, Eq. (J2H), without unitarity effects (denoted GRV94) using the 
GRV94 LO parameterization |23j in calculating the dipole cross section. The solid curves 
are the result including unitarity effects with the same GRV94 parameterization as initial 
input. The use of this parameterization is justified properly in Refs. 0, 124J. The dot-dashed 
curves are calculated with the dipole cross section, Eq. (J5TJ, using as input the GRV98 
parameterization for the gluon structure function. The aim of this comparison is to verify 
to what extent an updated parameterization can absorb unitarity effects. It is verified that 
at RHIC energy, the unitarity effects could be absorbed in the parameterization. However, 
at LHC energy the situation is quite different, and the results are completely distinct. The 
deviation is important mostly at large px, and as a general feature concerning the unitarity 
effects, those corrections are significant at large transverse momenta and are enhanced as 
the energy increases. 




01234 01234 

p T (GeV) p T (GeV) 



FIG. 4: The Drell-Yan dilepton transverse momentum pr distribution at RHIC (y/s = 500 GeV) 
and LHC (1/s = 14 TeV). The solid lines are the GM results including unitarity effects, the long- 
dashed ones are the curves using GRV94 for the gluon distribution (without unitarity effects) in 
the dipole cross section. The dot-dashed curves are the results obtained with the GRV98 gluon 
distribution (without unitarity effects) in the dipole cross section and the dotted ones are the 
results using the BGBK model. 
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FIG. 5: The Drell-Yan dilepton transverse momentum distribution at RHIC (y/s = 500 GeV) 
and LHC (yfs = 14 TeV). The curves have the same meaning as in Fig. @J At RHIC energy, only 
the GM and BGBK results are shown. 

As an additional comparison, we present curves from the improved saturation model, 
Eq. P4*|) (dotted lines in Fig. HJ): at RHIC energies and quite small p?, BGBK results 
overestimates the GM one; however, at high pt the BGBK model underestimates the GM 
predictions. At LHC energies, the BGBK underestimates the GM results. It is worth 
mentioning that until here the analysis has been performed for fixed values of mass and xp, 
which implies that the values of the variable x<i remain almost unchanged in the analyzed 
Pt interval. The unitarity effects studied here are calculated perturbatively, and thus they 
are more significant at small r. At small px, large r contributions are important and even 
dominate in that region, which does not allow to observe the saturation effects in a clear 
way. There, the confinement aspects of the process are more important. In contrast, at large 
Pt the main contribution comes from the small r region, which is sensible to the inclusion 
of unitarity corrections to the process. 

In order to perform estimates for more realistic values of the kinematical variables, one 
considers that the DY measurements at these colliders will be made predominantly in the 
central rapidity region 25], i.e. at xp — 0, instead of a very forward direction. In Fig. 
we present estimates for RHIC and LHC energies at xp = 0. For RHIC, the results for 
the BGBK model (dotted line) and for the GM approach (solid line) are shown, where the 
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deviations are larger at small pp. For xp = 0, the deviations due to the unitarity effects are 
smaller than for xp = 0.625, so only the GM distribution is shown, since the results with the 
GRV94 and GRV98 are almost the same as the one with the GM distribution. The results 
for LHC are also presented. The unitarity effects are smaller at xp = 0, because X2 is larger 
at midrapidity than in the forward direction. 

As a final investigation, the x^-integrated dilepton transverse momentum distribution is 
calculated and compared with the available data on pp reactions at y/s = 62 GeV and mass 
interval 5 < M < 8 GeV (CERN R209) Q] The results are presented in Fig. El with the 
solid curve denoting the Glauber-Mueller calculation, including the non-asymptotic valence 
content (GM + Reggeon), the dot-dashed line is the BGBK result [l| and the long-dashed 
line is the Glauber-Mueller calculation without non-asymptotic valence content (GM no 
Reggeon). The calculation using the improved saturation model shows only fair agreement 
with the experimental CERN data. Note however, that no reggeon part has been introduced 
for the BGBK model. In addition, the data shown in Fig. IHlwere integrated over all xp and 
therefore include contributions that are not taken into account by the dipole approach (see 
discussion in sect. EJ). The GM result, on the other hand, is in good agreement with 
the overall normalization and behavior presented by the data, when the non-asymptotic 
contribution is taken into account, even though no parameters have been adjusted to fit the 
data. The GM cross section overestimates the saturation model due to the inclusion of the 
non-asymptotic contribution. 



VI. CONCLUSIONS 



In this work, we investigated in detail the Drell-Yan transverse momentum distribution 
in the color dipole framework, and we analytically demonstrated that the dipole approach 
correctly reproduces partially the NLO parton model in the appropriate limit. In contrast 
to the cross section integrated over pp, the DY pp distribution opens a kinematical window 
where even large dipole configurations contribute. This can be verified by studying the 
weight function associated with the light cone wave functions for the process for different 
values of transverse momentum. Large partonic configurations have their maximal contri- 
bution at pp = 0. A remarkable feature of the dipole approach is the finite and well behaved 
property of the dilepton pp distribution at pp — > in a LO calculation. 
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2.0 3.0 

P T (GeV) 



FIG. 6: The Drell-Yan differential cross section on p^ at energy y/s = 62 GeV. The solid line is 
the GM result and the dot-dashed is the results using the BGBK model (without Reggeon part). 
The long-dashed line is the GM results without Reggeon contribution. 

The main motivation to pursue the dipole approach is that it provides a natural framework 
for the description of unitarity effects, that are not taken into account by the conventional 
parton model. Unitarity corrections are implemented in the dipole cross section, using 
the GM approach Q]. In addition, it is performed a comparison with the QCD improved 
saturation model of the dipole cross section 0. In g e„e r a,, the onitarity sections produce 
a reduction in the differential cross section, mostly at large pt- At LHC energies, the 
corrections are quite large and they cannot be reproduced by only using new adjustable 
parameterizations for the gluon distribution. 

In order to extrapolate the dipole approach to lower energies, a Reggeon contribution was 
introduced into the dipole cross section. This Reggeon part is proportional to the valence 
quark content of the target, meaning at high energies, i.e. RHIC and LHC, it is negli gibl e, 
although it is important in order to obtain a good description of the CERN ISR data 26] . 
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